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The connection between the out of equilibrium linear re- 
sponse function and static properties established by Franz, 
Mezard, Parisi and Peliti for slowly relaxing systems is ana- 
lyzed in the context of phase ordering processes. Separating 
the response in the bulk of domains from interface response, 
we find that in order for the connection to hold the interface 
contribution must be asymptotically negligible. How fast this 
happens depends on the competition between interface cur- 
vature and the perturbing external field in driving domain 
growth. This competition depends on space dimensionality 
and there exists a critical value d c = 3 below which the in- 
terface response becomes increasingly important eventually 
invalidating the connection between statics and dynamics as 
the limit d = 1 is reached. This mechanism is analyzed nu- 
merically for the Ising model with d ranging from 1 to 4 and 
analytically for a continuous spin model with arbitrary di- 
mensionality. 
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I. INTRODUCTION 

The off-equilibrium character of the time evolution of 
a system undergoing a phase ordering process, such as a 
ferromagnet quenched below the critical point, is clearly 
manifested by the aging property observed in the re- 
sponse function. If the system is cooled in zero field and 
left in the low temperature phase for a time t w before ap- 
plying an external field, for t w sufficiently large the time 
dependent magnetization displays a behavior of the type 



M(t, t w ) ~ M st {t -t w ) + M ag (t, t w ) 



(1) 



where M st (t — t w ) is a stationary time translation in- 
variant (TTI) contribution and the remaining term 
M ag (t,t w ) is the aging contribution obeying the scaling 
form 



M ag (t, t v 



t~ a M 



(2) 



A structure of the same type shows up also in the auto- 
correlation function given by 



C(t, ty 



Cstit t w ) + C ag I - — 



t 



(3) 



Behaviors like (Q) and @ are common features of slow 
relaxation and are the object of very intensive study es- 
pecially in glassy systems, with and without disorder [[!]. 

In the case of systems evolving via domain coarsening, 
structures of this type can be readily interpreted in terms 
of two independent variables responsible, respectively, of 
the fast thermal fluctuations within domains and of the 
slow out of equilibrium interface dynamics. The splitting 
of the order parameter into thermal and ordering com- 
ponents was used some time ago Q as the key ingredient 
in the theory of phase ordering. Therefore, the station- 
ary contributions in (|l]) and ([?]) are due to equilibrium 
thermal fluctuations in the bulk of domains, while the 
aging terms come from the remaining out of equilibrium 
fluctuations occurring at the passage of interfaces . 

In the study of glassy systems, along with the realiza- 
tion that in these systems the out of equilibrium proper- 
ties are of foremost importance, recently there has been 
a pair of developments which have further enhanced the 
interest in the dynamics of slow relaxation. The first has 
been the off-equilibrium generalization of the fluctuation 
dissipation theorem (FDT), first derived by Cugliandolo 
and Kurchan j^] in the context of mean field models for 
spin glasses. This amounts to the statement that for 
t w — > 00 the magnetization depends on the time vari- 
ables only through the autocorrelation function 



M(t,t w )=M[C(t,t w )] 



(4) 



and the deviation from the ordinary FDT can be ex- 
pressed through the so called fluctuation dissipation ratio 
(FDR) 



X(C) 



dM(C) 
dC 



(5) 



which obeys X(C) = 1 in equilibrium. The second is a 
theorem by Franz, Mezard, Parisi and Peliti (FMPP) || 
which allows to retrieve the structure of the equilibrium 
state from dynamic properties during relaxation. Under 
certain hypothesis, they have established the identity 



dX(C) 
dC 



= P(q) 



(6) 



C=q 



where P{q) is the overlap probability distribution in the 
equilibrium state M, This latter development is of par- 
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ticular significance, since it opens a way around the dif- 
ficulty of static computations for systems with complex 
equilibrium states. 

In this context, the phase ordering process in pure 
systems is of considerable interest as a simplified frame- 
work where the chain of connections aging-FDR-statics 
can be analyzed and tested. The main point is that 
the phenomenology of phase ordering displays the typ- 
ical features (|j) and (||) of slow relaxation, and that the 
structure of the equilibrium state is exactly known, thus 
allowing for a detailed investigation of the relation be- 
tween statics and dynamics. Work in this direction |||]]9) 
has led to the conclusion that the aging term in the re- 
sponse function does not play any role asymptotically, 
therefore relegating the connection between static and 
dynamic properties in the somewhat trivial bulk contri- 
bution. The argument is based on the statement that 
interface response comes only from the spins on the bor- 
der of growing domains, yielding the upper bound 



M ag {t,t w ) < pi{t w ) ~ L-\t w ) 



(7) 



1 / z 

where pi(t w ) is the interface density and L(t w ) ~ tj, is 
the typical domain size. This fits into the form (^) with 
a = l/z, where z is the growth exponent. 

However, this picture is at variance with exact analyt- 
ical results for the one dimensional Ising model 10 ll| in 
the limit of infinite ferromagnetic coupling (the reason for 
taking this limit rather than the zero temperature limit 
will be discussed in Section 6). In this case one finds 
the opposite situation, namely there is no bulk response, 
while the interface response obeys (0) with a = and 



M (—} = — arctan,/— - 1 (8) 

\t w J 7T V tw 

yielding a finite asymptotic value independent of t w 

lim M ag (t,t w ) = (9) 
Similarly, there is no stationary term in the autocorrela- 



tion function, while the aging term is given by |12 1 



(' '„,,, [ j- ) = - ai'csin 



(10) 



Hence, eliminating t/t w between (||) and ( [Toj ) one finds 
(Fig- 1) 



V2 

M (C) = arctan 

7T 



\/2cot (|c 



(11) 




FIG. 1. M(C) for the d = 1 Ising model with J = oo. 

showing that the response function obeys (Q) for any 
t w giving rise to a non trivial FDR which, however, leads 
to a violation of the connection (^) between static and 
dynamic properties. This result indicates that the ag- 
ing part of the response function for coarsening systems 
might not be as simple as (Q). In order to address this 
problem jlj] , we have analyzed the behavior of the in- 
terface contribution to the response function as dimen- 
sionality is varied, finding through simulations for Ising 
spins and a phenomenological model for continuous spins 
that the interface contribution is indeed less trivial than 
hitherto believed. We find that the scaling form (|J) holds 
with a scaling function and an exponent a which depend 
on dimensionality, providing a unified coherent picture 
for the diverse behaviors observed at different dimen- 
sionalities. More specifically, we find that d c = 3 is the 
critical dimensionality such that: 

i) for d > d c the response is actually due only to the 
polarization of the spins at the interfaces making (Q) to 
hold 

ii) for d < d c there is a new and non trivial behavior 
of the response function due to the competition in the 
motion of interfaces between the drive of the curvature, 
aiming to minimize surface tension, and the drive of the 
external field, aiming to minimize the magnetic energy 
of domains. 

The paper is organized as follows. In Section 2 gen- 
eral concepts about the structure of phase space and time 
evolution are reviewed. Section 3 and Section 4 are de- 
voted, respectively, to the relaxation process dominated 
by the fast degrees of freedom leading to equilibration 
and to the phase ordering process which is, conversely, 
dominated by the slow out of equilibrium degrees of free- 
dom. Section 5 contains a short account of the FMPP 
scheme for the connection between static and dynamic 
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properties. Section 6 and Section 7 contain results, re- 
spectively, for the Ising model in d = 1 and in higher 
dimensions. The model for continuous spins is presented 
in Section 8 and concluding remarks are made in Section 
9. 



II. STRUCTURE OF PHASE SPACE 

Let us consider a spin system with hamiltonian W[sj], 
for instance the ferromagnetic Ising model, in contact 
with a thermal reservoir at the temperature T. Below 
the critical temperature T c configuration space breaks 
up into ergodic components |]l5| 



n = (u a n a ) un b 



(12) 



where by fl a , with a = ± we have denoted the basins of 
attraction of pure states and by fib the boundary between 
them The Gibbs state 



Pg(w) = — ' 



(13) 



is the mixture of the two broken symmetry pure states 

Pg{u) = w + p+(uj) + w-p-{uj) (14) 

where Z = XLen e^^ n< - u \ to — [s,] is a spin configura- 
tion, the pure states are given by 



, if u> ^ Q a 



(15) 



and iu+ = W- — 1/2. In each pure state there is sponta- 
neous magnetization 



m± = 



N ^ 



Si)± = ±m T 



(16) 



and a finite correlation length £r which does not depend 
on the sign of the state and is related to the relaxation 
time within the pure state by 



A/z 
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(17) 



In what follows t will characterize the time scale of 
fast relaxation. 

Alternatively, defining the overlap of two configura- 
tions by 



(18) 



the structure of a state p(u>) can be characterized through 
the probability that Q(u),u>') takes the value q when 
lj and u>' are configurations of two independent copies of 
the system 



P(g)=^p(o;)p(a; / )5(Q(^^)-9)- 



(19) 



Using (|T^) and (|l6|), the overlap probability function in 
the Gibbs state is given by 

P G (q) = (w% + w 2 _)6(q - m\) + 2w+W-5(q + m\) (20) 

where the mixed character of the state is revealed by the 
presence of the second ^-function in the right hand side. 



For future reference, notice that from (19) follows 



dqP(q)q 



N 4^ 



(21) 



Let us now consider the instantaneous quench process, 
where the system is initially prepared in some initial state 
Poiijo) and, at the time t = 0, is put in contact with the 
thermal reservoir at the temperature T < T c . Taking 
t > t, the measure over is given by 



t) = po(Q a )p a (ui) + pb(ui, t) 



(22) 



where po(fl a ) = J2uj^a a /°o(w) and pb(u), t) is the measure 
over the boundary. Similarly, for the joint probability at 
times t > t' > t we may write 

p(uj't', Ut) = ^2 Po(Qa)p a (u' , L0,t- t') + pb{u' t' , Ujt) 

a. 

(23) 

where p a (u> ,u>,t — t') is the TTI pure state joint prob- 
ability. From (|2^) and ( |23| ) it is quite clear that the 
properties of the system following a quench below T c are 
sensitive [[l5| to the choice of the initial condition po(u>), 
specifically to the weight given at the time t = to the 
different components. 

At the level of the observables of interest, like magne- 
tization m(t) — (si(t)} and correlation function C{\ i — j \ 
,t,t') — (si(t)sj(t')} — (si(t))(sj(t')), where space trans- 
lation invariance is assumed to hold, the above results 
translate in the following way. From (|22| ) follows that for 

t > T 



m(t) 



po(tt a )m a + m b (t) 



(24) 



where m a is the equilibrium value of the magnetization 
in the pure states given by (|l6|). Next, assuming that on 
the boundary uib(t) — 0, for t > t' > t and from ( p3| ) we 
have 



C(\i-j\,t,t') = 



^2po{n a ) 



C ps (\i-j\,t-t') 



+ C b {\i-j \,t,t') + Am 



(25) 



where 
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C ps (\i-j \,t-t') = ( Sl (t) Sj (t')) c 



(26) 



is the TTI correlation function in the equilibrium pure 
states which, for pure states related by symmetry, is in- 
dependent of a. Furthermore 

C b (\i-j\,t,t') = (s i (t)s j (t')) i (27) 

is the correlation function on the boundary and 

-i 2 



Am = p (n a )m 2 a 



y^p (Sl a )m a 



(28) 



gives the fluctuation of the magnetization over pure 
states. Properties of the pure state correlation function 
which will be needed in the following are 



C ps (\ i-j\=0,t-t' = 0) = l-m 
where we have used sf = 1, and 

C ps {\i-j |,t-*')=o 

for | i - j | > £ T or t - t' > r. 



(29) 



(30) 



III. FAST PROCESS: RELAXATION TO 
EQUILIBRIUM 

Let us now adjust the initial condition of the quench 



in order to have relaxation to the Gibbs state (13). 
From (14) and (FL5|) the Gibbs state is the mixture of 



pure states with weights w a = Z a /Z . Hence, according 
to (^H) relaxation to the Gibbs state can take place only 
if the initial condition is such that 



Po(ttb) = 0. 



(31) 



With such an arrangement, equilibrium is reached in the 
time scale r. Having now ^ Q po(Q a ) = 1, with po(£l a ) = 
1/2 independent of a, ( |24| ) and (25) simplify to m(t) = 
and 

C(\i-j\,t-t') = C ps (\i-j\,t-t')+m 2 T (32) 

implying the no clustering property C(\i — j\,t — t') > 
m?p. This behavior is illustrated in Fig. || depicting the 
autocorrelation function C(t — t') = C(\i — j\ = 0,t — t') 
for the d = 2 Ising model quenched to T = 0.969T C with 
the initial condition po(o-i) — \8{ui — lo + ) + ^5(uj — w_) 
concentrating the measure on the bottom of the basins 
of attraction, where lj + = [s^ = 1 , Vi] and w_ = [s.; = 
-l,Vi]. 
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FIG. 2. Autocorrelation function for the d = 2 Ising model 
with J = 1 quenched to T = 0.969T C (m| = 0.61584) from 
the initial condition po{u>) = \5{lo — ^+) + \S{lu — lu-) with 
W4. = [si — l,Vi] and = [si — — l,Vi]. 

An observation important for what follows is that the 
form ([32J) of the correlation function corresponds to the 
splitting of the spin variable into the sum of two statis- 
tically independent components 



Si{t)=1pi(t) 



(33) 



where i/ji(t) obeys the statistics of equilibrium thermal 
fluctuations in the pure state with expectations 



(V»i(i))a = 



C ps (\i-j\,t-t'). 



(34) 



Instead, a is a time independent random variable, the or- 
dering component, which takes the values m a with prob- 
abilities p(m a ) = po(fl a ) determined by the initial con- 
dition. Denoting the latter average by an overbar we 
have 



-5 = 

\ <J 2 = 

and (j|^) can be rewritten as 

C(\i-j\,t-t') = (i H (t)i; j (t , )) c 



(35) 



(36) 



The next step is to study the response of the system 
to a perturbation. Suppose that at the time t w after the 
quench the hamiltonian Ti.(uj) is changed into 



H h (oj) = U{uj) - Hi(u) 



(37) 



where 



(38) 
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is an uncorrelated gaussian random field (RF) with ex- 
pectations 



E h {hi) = 
E h (hihj) = h^Sij. 



(39) 



In the RF Ising model the lower critical dimensionality 
is raised from d^ = 1 to d^ = 2. Hence, for d > 2 the 
component structure ([l2]) of configuration space is not 
modified by the presence of the RF. Due to the pres- 
ence of the perturbation the system will relax to a new 
equilibrium state 



Po(&a)pa,h(u) 



which is neither the perturbed 



PG,h( u ) 



(40) 



(41) 



nor the unperturbed Gibbs state (p^3[), since the pertur- 
bation is present in the pure states, but not in the weights 
po(Q a ). We will be interested in the linear response to 
the perturbation, since already in the simple context of 
fast relaxation it is possible to identify some of the basic 
elements of the connection between static and dynamic 
properties to be discussed in Section 5. Let us then con- 
sider the staggered magnetization defined by 



T 



M{t-t w ) = j^E h [{Hi{t)) h ] 



T 



E h 



. i 

(42) 



where (-)h denotes the thermal average for a given RF 
realization. Taking t w > r, i.e. switching on the pertur- 
bation after the unperturbed system has reached equi- 
librium, this quantity depends only on the time differ- 
ence. Since by definition (ipi(t)) = 0, we may write 



(si(t))h = Oi{t) where, due to the RF, the ordering 
component a acquires a site and time dependence. Ex- 
panding up to linear order in the field and recalling that 
the unperturbed average a vanishes we have 



= ^x{\i - - tv>)hj + 0{h 2 ) 



(43) 



where from ( |40| ) 

X(|*-j'|,*-*w) = Ys M^a)Xa(|« - j\,t - t w ) (44) 

OL 

and for pure states related by symmetry X a (|i — j\,t— t w ) 
is independent of a. Inserting (^) and ( |44| ) into ( f42| ) we 
obtain 



M(t-t w )=T Xa (t-t w ) 



(45) 



where Xa(t ~ t w ) = Xu(\i - j\ = 0,t - t w ). Hence, 
lim^oo M(t — t w ) = M a = T\ a where, using (p9|), the 
static susceptibility in either one of the pure states is 
given by 



1 



Xa 



[1 



ji "TJ 

This can also be rewritten as 



where 



M a = l- J dqP a (q)q 



P a (q) =S(q- rr? T ) 



(46) 



(47) 



(48) 



is the overlap probability function of the unperturbed 
pure states. We call the attention here on the point that 
this result is different from what one obtains computing 
the susceptibility from the perturbed Gibbs state dill), 
which differs from ([lo]) because the RF dependence enters 
also in the weights w a> h- In that case, in place of ( |i^ ) 
one finds 



XG 



(49) 



where (-)g denotes the average over the unperturbed 
Gibbs state. Recalling (p|), this can be rewritten as 



XG 



1 

T 



I- J dqP G (q)q 



(50) 



where now P G (q) is the overlap probability function (|2C 
of the Gibbs state. However, the form (49) or jscj) of the 
susceptibility cannot be reached dynamically. That is, 
by switching on the perturbation at the time t w after the 
quench, the t — > oo limit of the staggered magnetization 
is given by ( f46|) and not by ( fl9"| ) , which gives \G = 1 /T 
since (si)o = 0. 

Next, let us turn to the dynamical side of (Eq) and let 
us show, for pedagogical purposes, that ( (46| ) and ( fl7|) can 
also be obtained from a dynamical object like the linear 
response function 



R(t-t') = 



S(si(t)) 
Shiit') 



(51) 



h=0 



without resorting to knowledge of the equilibrium state. 
The response function Xa(t~t w ) entering (^) and R(t — 
t') are related by 



Xa{t - tw) 



J dt'R(t-t') 



(52) 



and, given that the unperturbed system is in equilibrium, 
the linear response function and the pure state autocor- 
relation function are related by the FDT 
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R(t - 1') 



1 dC ps (t-t') 



T 



dt' 



(53) 



Since the constant term in ( p2[) makes no contribution 
to the time derivative, we may replace C ps (t — t') by the 
full autocorrelation function C(t — t'), and inserting ( |53| ) 
in (52) we find 



M{t-t w ) = 



f 



dq=[l- C{t - t u 



(54) 



This shows that when FDT holds the time dependence 
of M , or is entirely absorbed in the linear dependence 
on the autocorrelation function (Fig. |J). 
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FIG. 3. M(t — t w ) against C(t — t w ) for the same quench 
as in Fig. |^. The horizontal dashed line represents the con- 
tinuation of M(C) into the unphysical region C < m^p. 

From ( |54] ) M(t — t w ) reaches the equilibrium value |f6| ) 
as C(t — t w ) reaches the lower bound rn^. Even though 
C(t — t w ) cannot fall below this value, we may extend the 
dependence of M on C into the unphysical region C < 
mip (horizontal dashed line in Fig. |3|) by rewriting ( [54] ) 
as 



M{C) = [ 9{q- m 2 T )dq 
Jc 



(55) 



where 9 denotes the Heaviside step function. Integrating 
by parts we find 



M(C) = [1 - C9(C - m%) 
which yields 



dqS(q — m^p)q (56) 



•W^Hm -I—' '^ti?- 1 (57) 



1 — rrifp , for C < mip. 



Taking C = 0, from @ we find M a = 1 - £ dqS{q - 
m^)q recovering pTj). In order to understand this result, 



notice that ( |55| ) c & n be regarded as obtained from (|5 
with the FDT in the modified form 



R(t - t') = 



X[C(t-t')} dC(t-f) 



T 



or 



(58) 



where X(C) = 9{C - m\) is the FDR introduced in (|). 
Rewriting (B9) as 



M(C) = [1 - CX(C)} - [ dq^^q 
Jc dq 



(59) 




and taking again C — we find M a = 1 — J Q dq dx J-^ q. 
Comparing with (^) then we find the relation (^) in the 
form 



dX(C) 
dC 



P a (q) 



(60) 



C=q 



showing that the piece of static information contained in 
P a {q) is encoded into the relaxation properties through 
the FDR. Although this may seem an artificial exercise, 
it will turn out useful in the understanding of the con- 
nection between static and dynamic properties in the less 
trivial context of slow relaxation. 



IV. SLOW RELAXATION: PHASE ORDERING 

In the previous Section we have analyzed a quench pro- 
*cess which yields equilibration in the Gibbs state within 
the time scale r. In order to achieve this the initial con- 
dition had to be chosen according to (|3l]). Now we turn 
to phase ordering [fL7f , where equilibrium is not reached 
within any finite time scale. We will find out that in or- 
der to have a phase ordering process the initial condition, 
in a sense, must be opposite to (|3l"| ) with 



Po(^a) = 
po(^h) = 1. 



(61) 



Nonetheless, the fast equilibration process of the previous 
Section will turn out to dominate the short time behavior 
of phase ordering. 

In order to assess how the phase ordering process fits in 
the general scheme of Section |J, we relay on the behavior 
of the correlation function. Typically, the initial state is 
taken as the infinite temperature equilibrium state 



(62) 



with p(si) = 1/2 which yields the uniform measure 
po(u>) = 2~ N . Taking the shortest time after the quench 
t' sufficiently larger than r, the observed behavior of the 
correlation function is well represented by the sum of two 
contributions 
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C(\i-j\,t,t') = C pa (\i-j\,t-t') + C ag (\i-j\,t,t') 

(63) 

where the first one is TTI and coincides with (^6|), while 
the second one displays aging through the scaling behav- 
ior E7p| 



C ag (\i-j\,t,t') 



m T F ag 



J\ 



L(t') 



L(t) \ 
L(t')J 



(64) 



The characteristic length L(t) grows with the power law 
L(t) ~ t x l z where z = 2 for non conserved order param- 
eter, as it will be considered in this paper. The scaling 
function F ag (x,y) has the properties 



F ag (x, 1) ~ e 



Fag(0,V) 



V 



x , for x > 1 
A , for y > 1 



(65) 



with A > 0. This shows that, contrary to the previ- 
ous case, now C(\i — j\, t, t') becomes smaller than m\ 
and eventually vanishes when \i — j\ or t — t' becomes 
large. This, in turn, implies that with the uniform initial 
state (^) condition (|6l| ) is realized, otherwise from (|25|) 
follows that it would not be possible for C(\i — j\,t,t') 
to vanish at large distances or large time separations. 
Therefore we must have 



C(\i-j\,t,t') = C b (\i-j\,t,t') 



(66) 



revealing that the structure ( |63| ) reflects a property of 
the evolution over the boundary flj,. In the simplest case 
of the ferromagnetic system with two pure states, as we 
are considering, it is well known that the time evolution 
of configurations is given by the coarsening of domains of 
the two opposite equilibrium phases. Taking t' S> r in or- 
der to separate the time scales of fast and slow dynamics, 
within each domain the system is in equilibrium in either 
one of the two pure states p±{u). Putting this together 
with (|6^), since in the time regime t — t' <C t' and for 
short distance C ag — m^, we have that over short time 
and short distances the correlation function is indistin- 
guishable from (]3^). Namely, we may regard the phase 
ordering process as a fast relaxation process of the type 
considered in the previous Section, followed by a quite 
different and much more slow relaxation taking place on 
the time scale set by t'. This suggests the separation of 
the spin variable into a fast and a slow component, by 
generalization of the split (^) in the form jj~ 



i(t) = 1pi(t) + (Ti(t) 



(67) 



where, again, ipi and Ui are two statistically independent 
variables. By analogy with (|33|), we define the slow 
ordering component by 



(Ji(t) = ±m T 



(68) 



according to the sign of the domain the site i belongs to 
at the time t. Then, the statistics of the one time prop- 
erties of at (t) is determined by the relative occurrence of 
domains of either sign, which yields the time independent 
probability p(o"i) = 1/2 and the expectations 



(Jiit) = 
a?(t)=rr, 



(69) 



as in (|35|). However, since crj(t) changes sign when- 
ever an interface passes through the site i, the two times 
statistics is determined by the out equilibrium interface 
motion. With this choice for the ordering component, ipi 
represents the fast thermal fluctuations in the bulk of do- 
mains with the statistics of the equilibrium pure states, 
which is independent of the sign of domains. From (|67~ 
then follows 



C(\i - j\,t,t w ) = (ipi(t)Tpj(tw))« + <Ti(t)<Tj(t w ) (70) 
and comparing with ( |63| ) we can make the identifications 
(ipi(i)^j(tw)) a = C ps (\i-j\,t-t w ) (71) 



CJi(t)aj(t w ) = C ag (\i - j\,t,t w ). 



(72) 



After surveying the unperturbed phase ordering pro- 
cess, let us go over to the behavior of the staggered mag- 
netization ( (42| ) after switching on the perturbation ( |3"8| ) 
at the time t w . As mentioned above, with the RF the 
lower critical dimensionality is — 2. We will assume 
that the external field is so small that the bound on the 
size of domains imposed by the Imry Ma length £(^o) 
for d < 2 is much larger than the size of domains L(t) 
in the time region of interest. Therefore, we shall deal 
with coarsening, irrespective of dimensionality Using the 
split ( |67| ) and following the argument of the previous Sec- 
tion we have again (si(t)}h — Oi{t) . Expanding up to 
linear order we generalize (^) by writing 

^dt) h = ^XB{\i - j\,t - t w )hj +xi(\i~ j\,t,t w )hj + 0(h 2 ) 

3 

(73) 

where the integrated response function has been sepa- 
rated into the sum of two contributions. The first one 
accounts for the change in the magnetization in the bulk 
of domains and, due to the separation of the time scales of 
fast and slow relaxation, is TTI. In other words, this con- 
tribution ignores the existence of interfaces and therefore 
under all respects is the same as the integrated response 
analyzed in the fast relaxation process of the previous 
Section, i.e. XB(t~t w ) — Xa(t — t w ). Instead, the second 
contribution accounts for the extra response due to the 
existence of interfaces and is not TTI. Inserting in (|4^), 
in place of (B3) we now have 
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M(t, t w ) = T XB (t - t w ) + T X i(t, t w ). 



(74) 



By definition, the bulk contribution obeys FDT and 
therefore, following the discussion of the previous Sec- 
tion, is related to the autocorrelation function by (57), 
with the difference that now the region C < is not 
unphysical. 

For what concerns the interface contribution, in the 
first time regime with t — t w <C t w interfaces can be re- 
garded as static and it is quite reasonable to take xi pro- 
portional to the interface density xi(t,tw) — Pi{t w ) ~ 
L (t w ). The question is what happens in the aging 
regime t — t w 3> t w dominated by interface motion. The 
usual argument (|^,|[ is that %i keeps on being propor- 
tional to the interface density 



Xi{t,tw) ~ Pi{t) 



(75) 



and therefore eventually disappears like i _1 (i). If so, 
it is clear that by taking t w large enough the interface 
contribution can be made negligible with respect to the 
bulk contribution, leaving ( p7|) to account for the rela- 
tion between the whole response and the autocorrelation 
function. However, as explained in the Introduction and 
as we shall see later on, the interface response function 
turns out to have more structure than what j7q ) allows 
for. In particular, there is a dependence on space dimen- 
sionality which cannot be accounted for by (T7q). 



V. STATICS FROM DYNAMICS 

Let us now come to the problem of detecting the struc- 
ture of the equilibrium state from the properties of the 
linear response function in the off-equilibrium regime. 
This requires, first of all, the generalization of FDT out 
of equilibrium. As we have seen in Section 3 when FDT 
holds the time dependence oi x(t — t w ) is absorbed into 
the dependence on C(t — t w ). In the FDT generaliza- 
tion proposed by Cugliandolo and Kurchan || , as stated 
in the Introduction, this holds also in the aging regime 
postulating that for large t w 

x{t,t w )=x[C(t,t w )}. (76) 

In order to establish the connection between the FDR ^) 
and the static properties, FMPP have considered the 
general case in which at the time t w the hamiltonian is 
changed into TL j [s] = TLq [s] — TL P [s] with a perturbation 
of the form H p [s] = J2 tl <..<i p Jh..i P s ii- a i p where the 
couplings Ji x ..i are independent gaussian variables. By 
considering the behavior of the expectation Ej[(H p (t)) j] 
they have derived the connection between the FDR and 
the overlap probability function of the equilibrium state. 
Here, for simplicity, we reproduce the main steps of the 
argument in the particular case of p = 1, when the per- 
turbation Hp reduces to the form (|3^) , referring to || for 
the treatment with arbitrary p. 



The expectation entering in the definition of the stag- 
gered magnetization can be written as 

Eh[(Hi(t)) h ]=E h J2p h {u ) t)H 1 {u) (77) 



where ph(u>,t) is the probability distribution evolving 
with the hamiltonian ([37]). Using the fact that hi are 
independent gaussian variables and integrating by parts 



E h [(H 1 (t)) h ]=h 2 E h 



8 



(78) 



The same quantity can be evaluated dynamically in 
the Martin-Siggia-Rose formalism Q and, assuming 
that (76) holds, one finds 



E h [(Hi(t))h] =fNE h 



1 — Ch{t, t w )Xh(Ch(t, t w )) 
(79) 



, dX h (g) 
dq — q 

C h (t,t w ) d Q 



where Xh and Ch are, respectively, the FDR and the 
autocorrelation function in the perturbed system. Tak- 
ing the t — ► oo limit and using lim^oo Ch(t, t w ) = 0, 
from @ and @ follows 



E h 



x ^ d 

4- dhl {s 



N 








1 - E h ( 


Jo 


~ f 





1 dX h (q) 
dq — q 



dq 



(80) 



where (-)h,oo denotes the average with the probability 
distribution 



lim ph(w,t) = ph(u, oc). 

t — ►OQ 



(81) 



Therefore, what we have up to this point is that un- 
der assumption (|76]) there exists a relation between the 
susceptibility in the state (^Tj) and the first moment of 
the FDR in the perturbed system. In the general case 
one has the same relation between the generalized sus- 
ceptibility with respect to Ji 1 ...i p and the p-th moment of 
the FDR. In order to go further on, more must be known 
about the properties of oo). In the context con- 

sidered by FMPP ph(y>-, oo) coincides with the perturbed 
Gibbs state (|4l|), from which follows 



d 

dhi 



>G.h 



1 

T 



>G,h 



(82) 



Inserting this in the left hand side of (|30|) and using ( |2l| ) 
one finds 



E h J d q P G , h (q)q = E h f dq^j^q. 



(83) 
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From this and from similar relations for the higher 
moments one may establish the identity PG,h(q) = 
f q X h {q) which yields P G (q) = ±X(q) where P G (q) = 

lim/j^o Pgji(q) and X(q) = lim^o X h (q). Eventually, 
after establishing under what conditions Po(q) and X(q) 
may be identified, respectively, with the overlap function 
Pg(q) of the unperturbed Gibbs state and with the FDR 
X(q) of the unperturbed dynamics, one has the connec- 
tion between the unperturbed statics and dynamics in 
the form 



P G (q) = 



dX(q) 
dq 



Here, we call the attention on the fact that to estab- 
lish (|8^) it is essential that ( |8^ ) holds in order to use (21). 
Furthermore, the derivative with respect to hi in the left 
hand side of (|8(]) acts according to how the RF enters in 
the asymptotic state (|3l| ) . Therefore, if instead of reach- 
ing the Gibbs state ( |41[ ) the asymptotic state Ph(u, oo) 
coincides with the state (|40|), in place of (p3) one finds 



Eh / dqP a . h (q)q = E h 



1 dX h (q) 

dq — q 

dq 



(85) 



where P a ^h(q) is the overlap function of the perturbed 
pure state. Hence, following through the argument il- 
lustrated above, in place of ( |S4|) one concludes that the 
unperturbed FDR is related to the overlap function of 
the pure unperturbed state 



P a {q) 



dX(q) 
dq 



(86) 



In summary, the static information contained in the 
FDR depends on how the perturbation H. p enters in the 
asymptotic state (|8l|). We must now see in what form 
the scheme applies to the phase ordering process. Sup- 
pose that the interface response function Xiifi t w ) can be 
neglected in ([[I]) for t w sufficiently large. Then, as ex- 
plained in the previous Section, assumption ( [76] ) is veri- 
fied and M(C) obeys (^7|) leading to ( |60| ) which coincides 
with d86|). This is confirmed by numerical simulations on 
the Ising model for d > 2 fip| which show evidence for 
convergence toward the form (571) in the parametric plot 
of M versus C as t w becomes large. Therefore, a behav- 
ior of the type (|37|) is the signature that the interface 
contribution to the response function is negligible and 
the phase ordering process behaves as the fast relaxation 



process of Section 111 



VI. ISING MODEL D = 1 

In this Section and the next one we analyze the lin- 
ear response in the off-equilibrium dynamics of the Ising 
model beginning from the one dimensional case where 
analytical results are available. 



The system is defined by the hamiltonian with nearest 
neighbor interaction H(uj) = —JJ2i s i s i+i: where J > 
is the ferromagnetic coupling. From equilibrium statisti- 
cal mechanics we know that the equilibrium correlation 

function behaves as (siSj)c = e 5 ^ and the correlation 
length is given by £t = — [lntanh (y)] . Therefore er- 
godicity is broken for K = J/T — oo, which requires 
either T = for J < oo or J = oo and T arbitrary. 
Solving dynamics with K = oo, the two time correlation 
function obeys the form (^53|) where only the aging contri- 
bution (|64|) is present with = 1. The explicit form of 



(84) the autocorrelation function is given by (10). The reason 



for the absence of the TTI contribution C ps (\i— j\, t— t') is 
clear since with K — oo the ordering component 0$ = ±1 
coincides with Sj and ipi in (^) vanishes identically. 
Namely, in the quench with K — oo domains are formed, 
phase space motion takes place over Of, , as demonstrated 
by the aging behavior of the correlation function, but 
thermal fluctuations are absent within domains leading 
to the absence of the TTI contribution in the autocorre- 
lation function. 

It is interesting to consider also what happens in the 
quench with K < oo. In this case there is no ergodic- 
ity breaking. Solving dynamics one finds the generalized 

scaling form @ C(\i - j\, t, if) =F^,±$ ) ,ig± 

with the limiting behaviors 



F\ 



m_ L(f) 

.£(f)'L(t')' 
p ( \±^i t-l_ 



L(t> 



for - > 1 

for — < 1 and - < 1 



(87) 



where the equilibration time r is defined by ([17]) with 
z = 2. The meaning of ( |87| ) is, first of all, that after 
a finite time t equilibrium is reached. Hence, taking 
t'/r > 1 one observes the TTI behavior pertaining to the 
stationary dynamics in the equilibrium state. However, if 
t although finite is sufficiently large to allow for t/r <C 1, 
then in the time regime (t— if) <^ t one observes the same 
behavior as in the K — oo case. This can be understood 
considering that for K large t ~ e AK is the characteristic 
time needed to overturn one spin originally aligned with 
both of its neighbors. Then, taking if <S t means that 
one starts to look into the system when domains are still 
much smaller than Immediately after and as long 
as t < t growth continues as in the K = oo case, i.e. 
without thermal fluctuations within domains. Thermal 
fluctuations do come into play only for t > r. When this 
happens, the creation of defects by thermal fluctuations 
balances the losses due to interface annihilation and leads 
to a halt in domain growth and to the establishment of 
thermal equilibrium. 

Let us see what happens upon applying the RF at the 
time t w > 0. Recall that the master equation for the 
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system evolving with the Glauber spin flip dynamics is 
given by 



dp (w,t) 
dt 



= J2 M-Si)p (RitJ, t) - w( Si )p (w, t)] (88) 



where Rita is the configuration lu with the i-th spin re- 
versed and uu(si) is the transition rate from lj to RiU) 
given by 



Hf xt Si ) 



where H\ nt ■ 
7 = tanh(2A 



Hf 



tanh(K 



(89) 



and 



2\^i— 1 ~t~ *^'t • 1 ) : 1 ' • • \ 

Taking A" = oo, let us first consider the 
behavior of ( |89| ) before switching on the external field, in 
the interval < t < t w . Since 7 = 1, we have i/| n< = 1 
if Si~is i+1 = 1 and H" lt — if Sj_is i+ i = —1. In the 
latter case the spin Si is at the interface between two 
domains of opposite sign with probability 1/2 per unit 
time to flip or not to flip. Conversely, in the former case 
the spin flips with probability 1 if it is not aligned with 
its neighbors, while it does not flip with probability 1 in 
the opposite case, when it belongs to the bulk of a do- 
main. As a consequence the only dynamics in the system 
is the unbiased random walk of interfaces, leading to the 
growth law L(t) ~ t 1 / 2 . 

After switching on the RF we are interested in the be- 
havior of the staggered magnetization (|42]), which is now 
convenient to regard as the correlation function between 
the external field and magnetization. Right at the start 
M(t w ,t w ) — 0, since at t — t w the RF and configura- 
tions are uncorrelated. However, for t > t w the transi- 
tion rate ( |S9| ) is modified by the factor involving Hf xt 
which introduces a bias in the flips at the interfaces in 
favor of the local field, while bulk flips remain forbidden. 
Accordingly, M(t,t w ) grows positive revealing that spin 
configurations tend to correlate with the field. However, 
a substantial difference arises in the two ways to produce 
the limit K = 00. If J < 00 and T = 0, M(t,t w ) rises 
from zero to a plateau value M (Fig. within a micro- 
scopic time to and M depends on t w according to (inset 
of Fig. § 



M(t w ) ~L-\t w ). 



(90) 
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FIG. 4. M(t,t w ) for the d = 1 Ising model with 
J = 1 quenched to T = for different waiting times 



(10, 10^,3.16-10' 



, 10 3 , 3.162 10 3 ,10 4 , 



3.162- 10 4 ,10 J from top 



to bottom). In the inset, the plateau value M(t m ) is plotted 

— 1/2 

against tw The solid line is the t m behavior. 

The reason for this is that after switching on the field 
the motion of interfaces continues for a short time until 
pinning takes place when the two opposite spins mak- 
ing up the defect at the interface fall into a defect of 
the same sign in the field configuration. At that point, 
since Hf xt = sign hi, the second factor in the right hand 
side of J89|) vanishes and the interface does not move 
anymore. This explains reaching the plateau and the de- 
pendence ( |90| ) of M on t w , since the contribution to the 
staggered magnetization comes only from the spins at 
the interfaces and goes like the interface density. Fur- 
thermore, since T = implies ho/T — 00 there is no way 
to access the linear response regime, no matter how small 
the external field is chosen. Conversely, if the K = 00 
limit is obtained taking J = 00 and with no restrictions 
on the temperature, while the unperturbed dynamics re- 
mains the same, interesting behavior is obtained with RF 
since i) T > allows to overcome pinning of the inter- 
faces letting also the bulk of domains to participate in 
the correlation of spin configurations with the RF and 
ii) the condition ho/T <^ 1 can be realized making ac- 
cessible the linear response regime. 

In the following we will concentrate in the linear re- 
sponse regime with J — 00 and ho/T <C 1. In this case 
the staggered magnetization is given by (||) . The first ob- 



servation, comparing with the general form (74), is that 
the TTI bulk contribution is missing, as expected from 
the above discussion on the absence of thermal fluctu- 
ations when K = 00. Hence, the result (^) is entirely 
due to the second contribution in ( |74|) , which however is 
totally different from the behavior ( |75| ) which one would 
expect on the basis of a straightforward interface con- 
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tribution. Rather than decreasing at large time, here 
M(t, t w ) displays a correlation of the spin configurations 
with the RF which grows with time, until reaching the 
finite limit (^) as t — > oo (Fig. Q). Having excluded a cor- 
relation effect due to thermal fluctuations or to spin po- 
larization at the interfaces, the increase in the correlation 
can only be due to the fact that interface motion is driven 
by the field. As we shall now see, the field driven mech- 
anism, without modifying the growth law L(t) ~ f 1 / 2 , 
induces large scale domain drift in order to optimize the 
position of the bulk of domains with respect to the RF 
configuration. It must be stressed that although involv- 
ing the bulk of domains, this contribution to the stag- 
gered magnetization has nothing to do with the bulk re- 
sponse function coming from thermal fluctuations, which 
are now absent. 

An insight on how the field driven mechanism works 
comes from the behavior of (||) for t — t w <C t w from 
which we find 



Xi (t,tw) = 



1 



2(t - ty 



(91) 



Since in this time regime the system can be regarded 
as a set of non interacting interfaces and the density of 
interfaces pi{t) at the time t ~ t w is given by L -1 ^), 
we can rewrite (^lj) in the form 



where 



Xl(t,t w ) = pl(t)Xeff(t,U 



Xeff(t,t w ) ~ (t - t w y- 



(92) 



(93) 



is the effective response associated to a single interface. 
Looking next to the large time behavior for t — t w t w , 
from (§) follows T X i(t,t w ) = l/y/2 - 0(tlt w )- x l 2 which 
implies that the effective single interface response follows 
the behavior (^3|) also for large time. 

In order to check on this interpretation we have com- 
puted analytically the behavior of Xi(t>tw) when in the 
system there is a single interface. This is done by prepar- 
ing the system at t = in a configuration lu containing 
only one interface at the origin, for instance taking Sj = 1 
for i < and s, = — 1 for i > 0. The computation of the 
response function can be carried out exactly (see Ap- 
pendix A) yielding 



X 



sing 



(94) 



which substantiates the above analysis. This unexpected 
result makes it clear that the interface response is not 
simply due to the polarization of the paramagnetic in- 
terfacial spins, but is a much more complex effect. At a 
generic time t the interface has explored a region of order 
(t — tw) 1 ^ 2 and energy can only be released by reducing 
the contribution to Hi (lu) coming from the visited region. 



This can be achieved if the interface motion produces a 
large scale optimization of the position of domains with 
respect to the random field. 

In summary, from the analysis of the linear response 
function in the quench of the d — 1 Ising model with K = 
oo, we have uncovered a new and non trivial behavior 



different from the pattern outlined in Section IV, which 



is the one usually expected for coarsening systems. The 
role of the bulk and interface terms in (|74|) is reversed, 
the response being dominated by the latter one with all 
the new features illustrated above. 



VII. ISING MODEL D > 1 

As stated in Section 4, for the Ising model in two and 
three dimensions there is numerical evidence that as t w 
becomes large the staggered magnetization displays the 
structure (|57| ). For this to happen, the interface contri- 
bution in (|74| ) must vanish and only the bulk contribu- 
tion must be left over. However, the exact analysis of the 
previous Section in the one dimensional case is a serious 
warning that the interface contribution might not disap- 
pear so easily as the argument ( f75| ) could make believe. 
Therefore, a careful analysis of the interface contribution 
is needed to find out whether the field driven mechanism 
of interface motion is at work also in higher space di- 
mensions. In order to do so it is necessary to give an 
anambiguous definition of which degrees of feedom must 
be cosidered interfacial. In particular, it is necessary to 
make clear whether the flip of a spin in the interior of 
a domain belongs to a bulk fluctuation or makes a new 
interface. The definition we adopt is the following. Con- 
sider two configurations lui+b and luj evolving from the 
same initial condition with the same thermal history, un- 
der the influence of the same external field, if present. 
While LUj + b evolves with the usual Glauber dynamics, lui 
is subjected to the restriction that flips of bulk spins are 
forbidden. A bulk spin is defined as being aligned with all 
its nearest neighbours. Then, all the spins surrounding 
topological defects in uii are considered interfacial spins. 
On the other hand luj + b contains defects which can be 
either associated to interfaces or to bulk fluctuations, the 
latter being determined by comparison with u>j. On the 
basis of this definition we have measured the interface 
response function Xi{tjt w ) by simulating the evolution 
of the ferromagnetic Ising model with nearest neighbor 
interaction for d = 1, 2, 3, 4 without bulk flips and start- 
ing from the high temperature uncorrelated initial condi- 
tion ([32|). We have included the simulation of the d = 1 
case, for which the exact analytical results of the previ- 
ous Section are available, in order to have a check on the 
numerical procedure. 

For the effective response associated to a single inter- 
face Xeff(t, tw) defined by (p2[), we have obtained (Fig. ||) 
the asymptotic behavior 
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Xeff(t,t w )~(t-t w ) a (95) 
where the numerical values of a are consistent with 



3-d 
4 





for d < 3 
for d > 3. 



(96) 



For d = 3 a power law fit yields a = 0.03. A fit of 
the same quality is obtained with the logarithmic form 
T Xe ff(t, t w ) = 0.33 + 0.066 • log(t - t w ). 
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FIG. 5. Xeff(t,tw) in the Ising model without spin flips in 
the bulk. For d — 1, 2, 3, 4, the temperature, waiting time and 
linear system size of the simulations are T = 0.48, 2.2, 3.3, 4.4, 
t w = 10 3 ,10 3 ,10 2 ,10 and L = 10 6 , 512, 200, 42 with 
J — 1 and averages over 170,6045,114,922 realizations. 
The dashed lines are power laws with the corresponding 
exponent a. For d = 3 the curve is well fitted by 
TXeff(t, t w ) = 0.33 + 0.066 ■ log(t - t w ). 

In order to check to what extent Xeff{t,t w ) is related 
to a single interface response, we have performed another 
set of simulations without flips in the bulk, starting with 
an initial condition containing one straight spanning in- 
terface in the middle of the system. The results of simu- 
lations are shown in Fig. ^ and indeed the data reproduce 
quite well the behavior of Xeff(t, t w ), except for d = 3. In 
this case the logarithmic behavior found for Xe//(^5 ^tu) is 
followed up to a certain time, beyond which the response 
speeds up considerably. The analysis of this behavior, 
for which we do not have an adequate explanation, re- 
quires numerical investigation at much larger times and 
goes beyond the scope of the present work. 




FIG. 6. The single interface response, obtained from simu- 
lations without flips in the bulk and with an initial condition 
containing one single flat interface. For d = 1, 2, 3, 4 tempera- 
tures of the simulations are T = 0.48, 2.2, 3.3, 4.4, with J = 1, 
t w = and averages over 151600, 19810, 69891, 3537 realiza- 
tions. The dashed lines are power laws with the corresponding 
exponent a. 

In summary, comparing Fig. |^ with Fig. ^, the identi- 
fication of Xeff(t,t w ) with the response associated to a 
single interface Xsing(t,t w ) is on the whole well founded. 
We must now extract the meaning of the overall behavior 
as dimensionality is varied. The main feature is that the 
power law growth ( |95| ) weakens as d rises from d = 1 to 
d = 3 and disappears above d — 3. The explanation of 
this behavior can be conjectured recalling that in the un- 
perturbed system interfaces perform an unbiased random 
walk in d = 1, while are curvature driven for d > 2. In the 
perturbed system in d = 1, as we have seen in the previ- 
ous Section, interfaces are field driven. This mechanism 
operates also in higher dimensions, except that it enters 
in competition with the curvature mechanism. The effect 
of the curvature becomes comparatively more important 
as d increases due to the increasing coordination number. 
According to our simulations d c = 3 is the critical value of 
the dimensionality, such that for d > d c the field driven 
mechanism is ineffective and interface motion is domi- 
nated by the curvature yielding a = 0. Therefore, for 
d > d c , interfaces respond only through the polarization 
of the interfacial spins. This response saturates to the 
asymptotic value over a microscopic time scale (curves 
for d = 4 in Figj5| and Fig.^). Conversely, for d < d Cl 
the field driven mechanism competes with the curvature 
yielding a > and this competition gets more effective 
as dimensionality is lowered. Finally, as the limit d = 1 
is reached from above, the curvature mechanism disap- 
pears and the effectiveness of the field driven mechanism 
becomes complete yielding a = 1/2. 
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VIII. CONTINUOUS SPIN MODEL 

The discussion of the previous Section makes clear that 
dimensionality plays a crucial role in determining the 
relative importance of the bulk and interface response. 
In order to clarify further this point, in this Section we 
present an analytical calculation of x/(f, t w ) in the frame- 
work of continuous spins which allows to vary dimension- 
ality freely. The approximations involved in what follows 
are too crude for an accurate quantitative reproduction 
of the results of the simulations. Nonetheless, even in this 
form, the continuous model is quite useful to capture the 
overall qualitative picture. 

The treatment of phase ordering with a continuous, 
scalar and non conserved order parameter <p{x, t) is usu- 
ally based |l7j on the time dependent Ginzburg-Landau 
equation 



= V J 4> + r</> 



+ r)(x,t) 



(97) 



where r and g are positive constants and rj(x, t) is a gaus- 
sian white noise with expectations 



(r](x',t')r](x,t)) = 2T6(t-t')5(x-x"). 



(98) 



The infinite temperature initial condition is imposed tak- 
ing 4>o(x) = 4>{x,t — 0) as an additional source of noise 
gaussianly distributed with expectations 



MZ) = o 

M%)Mx') = AS(x-x'). 



(99) 



In past work on phase ordering kinetics most of the effort 
has been devoted to the study of the scaling properties of 
the equal time correlation function, that is, in the present 
language, to the study of C ag (\ i — j \,t = t') in ( pjj ) 
and (|64|). For this purpose, thermal fluctuations are 



usually neglected eliminating the thermal noise in (97) 
Then, one deals with the equation 



dcp_ 
~dt 



V 2 + r<f) - 



(100) 



where the only source of noise is the initial condition 
(j>o(x). 

From the analytical point of view, one of the most 
successful tools of investigation of this problem has been 
the gaussian auxiliary field (GAF) approximation which 
goes back to the pioneering work of Ohta, Jasnow and 
Kawasaki jn|[l7| ■ The method is suited to study the late 
stage, after local equilibriu m w ithin domains has been 
achieved. In the case of Eq. (10C) this means that locally 
the order parameter must sit at the bottom of either one 
of the two degenerate minima of the potential satisfying 



4>{x, t) = ±mo 



(101) 



where mo = \fvjg is the T = equilibrium value of the 
order parameter. The idea at the basis of the GAF ap- 
proximation is that ( |101| ) can be implemented through 
a non linear transformation on an auxiliary field u(x, t) 
over which perturbative methods can be applied. Differ- 
ent versions of the approximation correspond to different 
realizations of the non linear transformation. Here we 
take the transformation of the Kawasaki, Yalabik and 
Gunton p3] type 



(x,t) = 



u(x, t) 



, u 2 {xt) 



(102) 



Then, if the non linearity of u(x*,t) is mild, un- 
bounded growth is allowed eventually yielding 4>{x, t) — 
mosign[u(x, t)] which enforces the requirement (101). In 
order to actually carry out computations, one has to solve 
the dynamics of u(x, t) induced by (100) via ( |102| ), as we 
shall do below. 

After this brief survey of the GAF method, let us go 
back to the equation of motion (^) including thermal 
fluctuations. A systematic treatment of this problem 
based on the Martin-Siggia-Rose formalism and on the 
split of the order parameter into ordering and fluctuat- 
ing components was worked out in Ref. Here, we 
follow the same idea working directly with the equation 
of motion. Let us split the order parameter as in (K37 



(x, t) = ip(x, t) + a(x, t) 



(103) 



with the aim of separating the fast thermal fluctua- 
tions from the slow ordering component. Inserting ( |103| ) 
in (1371) we find 



dtp da 
~dt + ~dt 



V 2 i(j + V 2 cr + rip + ra - gtp 3 
3gip 2 cr — Sgtfja 2 — 3ga 3 + T] 



(104) 



and let us decouple ip from a replacing the mixed terms 
by 3gip 2 cr — > 3g(ip 2 )a and igipa 2 — > 3gipa 2 . Further- 
more, let us assume that T is sufficiently lower than T c 
to justify the self-consistent linearization ip 3 — > 3(i/j 2 )ip- 
Stipulating that ip is driven by the thermal noise and 
that a is driven by the noise in the initial condition, wc 
obtain the pair of equations 



dip 
~dt 



W 2 ip 



r-3g(iP 2 (x,t)) -3ga 2 (x,t) 



and 



da 
~dt 



= \7 2 a +[r- 3.g(</> 2 (x, i)>] a - ga 3 



1p + 7] 

(105) 
(106) 



with the initial conditions ip(x,t = 0) = and a[x, t = 
0) = 4>q(x). Let us make the assumption, to be verified 
a posteriori, that ip is the fast variable with relaxation 
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time r. Defining r e 



3g(tp 2 ) eq and making the 



additional assumption that within the same time scale 
a(x,t) reaches local equilibrium with a 2 {x,t) = rrv^ = 



for t > t in place of (105) and ( |106| ) we may write 



^ = VV-&V + »7(2,*) (107) 

0(7 

— = \7 2 a + r eq (T-ga 3 (108) 
where the equilibrium correlation length £t is given by 



6 



2r e 



(109) 



From ( 103 ) follows that the autocorrelation function is 
given by the sum of two contributions as in (|6^) C(t, t') = 
C ps {t - t') + C ag (t, t') with the TTI piece 

— t ~ t ' 

C ps {t-t') = {i}{x,t)i){x,t')) = {ijj 2 {x)) eq e ^ (110) 
and the aging contribution 

C ag (t,t') = a(x,t)a(x,t'). (Ill) 



The latter one can be computed from ( ]10q ) using the 
GAF approximation with the non linear transformation 
of the type ( |102[ ) in which TOq is replaced by 



<j(x, t) 



u(x, t) 



1 



(112) 



From ( |110| ) indeed follows that ip(x,t) describes the fast 
equilibrating thermal fluctuations with the characteristic 
time t = 

Consider, next, the effect on the ordering component 
of an RF with expectations analogue to (53) 



E h [h{x)]=0 

E h [h(x)h(x')} = hl8(S-x'). 



(113) 



For t > t w the equation of motion (|108|) is modified into 



da 2 o , , » 

— = V cr + r eq a - ga + h{x) 



(114) 



while (107) for i/j(x,t) remains unaltered. In order to 
generalize the GAF approximation, notice that the ex- 
ternal field affects the transformation ( |ll2| ) in two ways: 
i) through the auxiliary field u(x, t) and ii) by shifting 
the saturation value imj- of domains. Therefore, we 
separate a bulk and an interface term writing a(x, t) — 
<jb (x, t) + a j (x, t) where 



a B (x,t) = / dx" xb{x - x', t - t w )h(x') 



(115) 



and 



ai(x,t) 



u h ( x, t) 



(116) 



1 



The latter one is constructed to account only for the effect 
of the external field on the interface motion, by keeping 
the saturation value at the unperturbed level imj, while 
the former takes care of the remaining perturbation on 
the bulk of domains. Hence, for the staggered magneti- 
zation the decomposition ( |74|) a pplies where, according 
to the discussion of Section IV, xb obeys FDT and is 
therefore related to the autocorrelation function by ((5?" 
Here we are interested in the interface contribution 



Xi(t,t v 



1 

n 



<7i(x, t)h(x) 



(117) 



and in order to compute this quantity let us go back 
to (114). Since we want to extract the dependence of 



Uh{x,t) on the RF up to first order, after substituting 
a = as + 07 and keeping into account that ub is a first 
order quantity we obtain 



dai 
~~dt 



V 2 cr/ 



eg CJ - .90/ + h(x) 



(118) 



where the effect of ob goes into a redefinition of the vari- 
ance ft.Q of the RF, which will be neglected in the fol- 



lowing. Substituting (116) for ai(x, t) and dropping the 



subscripts h, the equation of motion for the auxiliary field 
is given by 



du ~ 



• eq 1 - 



where 



and 



a'(u) 



a»{u) = -3-^ 



,2 1 



(119) 



(120) 



(121) 



Performing, next, a mean field approximation by keep- 
ing only the lowest order non linear contribution in each 
term and linearizing self-consistently, after Fourier trans- 
forming over space we find 



du(k, t) 

at 



-[k 2 + D(t)]u(k,t) + h(k) 



where 



D(t) = -r e? + ^(Vii) 2 . 



nr. 



(122) 



(123) 



Defining the linear response function by 
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(124) 



with 



y(t,0) = efo dsD(s) (125) 

the formal solution of ( |122| ) reads 

u(fc, t) = R(k, t, Q)u(k, 0) + Xu(k, t, t w )h(k) (126) 

where Xu{k,t,t w ) — J dt'R(k,t,t') is the integrated re- 
sponse function of the auxiliary field. Carrying out the 
self-consistent computation of D(t) (Appendix B), the 
large time behavior of Y(t, 0) is given by Y[t, 0) = 

d + 2 

a (t + jh 1 ) 4 wnere A is a momentum cutoff and a = 



3Ad 



Inserting in (124). from Xu(t,t w ) 

4r e ,m|,(83r)3 

(2ir)~ d J dkxu{k,t,t w ) and t w > 1/A 2 follows 

/ n . . d d + 2 f . . d + 2 ( 1 

Xu(t, t w ) = (4tt)- 2 1— J dt't'-— f t _ t' + — 

(127) 

Similarly, for the unperturbed autocorrelation function 
of the u field we find (Appendix B) 



ITT A[(t+^)(^ + ^)]^ 



«(£, *)«(£, tf) = — ^ IA ' . (128) 

a 2 (4^)<V 2 (t + <' + ^) d/2 



Now, in order to compute ( |117| ) w e m ake a further mean 
field approximation by replacing (|l!6j ) with 



<7/(x, f) = m,T 



Uhjx, t) 



(129) 



where S{t) = u 2 {x,t) is an unperturbed average. This 
gives 



Xi(t,t w ) 



and computing S(t) from ( |128| ) we get 

1 

2A 1 



5(*) = 6(*+^7a 



with b = 4re g" T and 



Xl (t,t w )=At^F( t -f, t -f 



where 



t w to 
T'7 



(130) 



(131) 



(132) 



(133) 



Here, to 

3d 



A 2 is a microscopic time and A 



4r cg (4ir) d 



Next, we may use the form (129) of the transforma- 
tion to compute also the aging contribution (111) to the 
autocorrelation function obtaining (Appendix B) 



C ag {t/t w ) = m\ 



d/4 



2t 



-d/2 



(134) 



For d — 1 the time ratio t w /t can be eliminated be- 
tween (132) and (134) yielding a parametric plot (Fig. [?]) 
of the response function versus C qualitatively similar to 
the one of the d = 1 Ising model in Fig. [|. In particular, 
in the large time limit we find the counterpart of (^|) 



lim Xi(t,t% 

t — >oo 



AF(0,0) 



2r. 



r(i/4) 



2tt 



(135) 




0.4 0.6 

c ag (t/t w ) 

FIG. 7. Parametric plot of \i against C ag in the continuous 
spin model. 

The interesting point now is to extract the behavior 
of the effective response function Xeff associated to the 
single interface and defined by (|92]). In the short time 
regime t — t w <C t w from (132) follows 



Xeff(t~t w ) = 



2A 



d-2 



2-d 



t ty 

to 



(136) 



which for to <C t — t w <C t w yields 

( 2A d ~ 2 



Xeff(t-t w ) = < 



d-2 



log(^) 

2A d ~ 2 ( t-t w \ ' 
t 2-d \ t J 



for d > 2 
for d = 2 



for d < 2. 



(137) 
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A similar behavior is obtained also in the large time 
regime t — t w > t w 



2A a 
d-2 



■A 



, for d > 2 

Xeff(t -t w ) = l 4A log , ford = 2 (138) 

I AF(0,0)t^ , for d< 2. 

Therefore, apart from a change in the prefactor taking 
place about t — t w ~ from (137) and (138) follows 
that both for short and large time Xeff obeys a power 
law as in (Pol) where, however, now 



a 



2-d 
2 





for d < 2 
for d > 2 



(139) 



and there is logarithmic growth for d = 2. The full time 
dependence of Xeff jt, t w ) obtained from the numerical 
computation of fll32| ) for different values of d is displayed 
in Fig. |^. Comparing Fig. ^| and Fig. |8|, the common fea- 
tures may be summarized stating that in both cases Xef f 
obeys the power law (^H|) and that there exists a critical 
value of the dimensionality d c such that the exponent a 
is zero for d > d c with logarithmic growth at d — d c . For 
d < d c the exponent a grows positive with decreasing 
dimensionality reaching the final value a = 1/2 at d = 1. 
The meaning of the critical dimensionality in relation to 
the growth mechanism has been discussed in the previous 
Section. The difference with the case of Ising spins is that 
now we have d c = 2 in place of d c — 3. This tells that, 
although qualitatively correct, the mean field approxi- 
mation developed above is not accurate enough to ac- 
count quantitatively for the competition between the field 
driven and curvature driven growth mechanisms. For in- 
stance, we find a domain growth law L(t) ~ t 1 / 2 also 
in d = 1, while the one dimensional continuous model is 
known to have logarithmic growth law plj] . Despite these 
shortcomings, the model reproduces the gross features of 
the response function as dimensionality is varied. 



10 



10 



10 



10 



10 



10 



FIG. 8. Xeff(t,t w ) in the continuous spin model with 
t w = 10 s . The dashed lines are power laws with the cor- 
responding exponent a. 



IX. CONCLUSIONS 

In this paper we have studied the behavior of the re- 
sponse function in non disordered coarsening systems un- 
der variation of space dimensionality. The results ob- 
tained are instructive on the applicability of the FMPP 
theorem in general. In order to clarify this point, let us 
go back to the form ([!]) of the staggered magnetization 







1/2 






// 
// 
// 
// 




d=1 / 


1/4 j 


1/2 


d=3/2^ 








d=3 





10" 



10 



10' 

t-t. 



10 



10 



M(C,t w ) = M st (C) + t- a M(C) 



(140) 



where we have used (g). As we have emphasized, the 
above pattern in the response function reveals the exis- 
tence of slow and fast degrees of freedom with widely sep- 
arated time scales. When P(q) is extracted from (14C) 
two basically different cases must be distinguished. If 
a ^ 0, as it is the case for coarsening systems with d > 1, 
the slow degrees of freedom for large t w make a negligi- 
ble contribution and the relevant information comes only 
from M st (C) yielding P st (q) = S(q - qsA), where q E A 



is the Edwards- Anderson order parameter (qEA 



0- 



Conversely, if a = one obtains an additional contribu- 
tion due to the slow degrees of freedom 



P(q) = P st (q) + P ag (q) 



(141) 



where P ag {q) = -d 2 M{C)/dC 2 \ c = q - This non trivial 
contribution appears in glassy systems and reproduces 
the expected pattern of replica symmetry breaking of the 
equilibrium state ||. However, this quantity appears also 
in the d = 1 Ising model with (Fig. ||) 



d 2 M{C) 



dC 2 



C=q 



[2- 



sm(f )] 2 



(142) 



and in this case it is not related to the equilibrium state. 
The general question then is: if during some relaxation 
process the response function takes the form ( 14C ) with 
a = 0, under what conditions M(C) contains informa- 
tion on the equilibrium state. The preceding analysis 
for coarsening systems suggests that the answer has to 
do with one of the hypothesis in the theorem, which re- 
quires the system eventually to equilibrate, and with the 
mechanism of slow relaxation. In glassy systems the time 
evolution proceeds toward equilibrium through decays of 
metastable states ||. This may take very long, but even- 
tually all degrees of freedom, including the slow ones, 
will equilibrate. The case of coarsening, instead, is qual- 
itatively different. Slow relaxation is not due to decay 
, of metastable states, there are no activated processes. 
Rather, there is a smooth reduction of defect energy, as 
motion in phase space takes place over the border, where 



1G 



the slow degrees of freedom do reduce in number but 
never equilibrate. Hence, in this case, M{C) is a prop- 
erty of an intrinsically out of equilibrium dynamics with 
no connection to any property of the equilibrium state. 



FIG. 9. 



d*M(C) 



for the d = 1 Ising model with J — oo. 



90 

dt 



= S7 2 (J) + T) 



(143) 



d(f>(k,t) 
dt 



-k 2 (j>(k,t) +T](k,t) 



(144) 



A similar behavior is displayed by the equal time corre- 
lation function C(t,t) = (<fi 2 (x,t)) for large time 



T(4tt) 



-A 2A a 



d-2 




C(t,t) 



T(4 7r )-# log (|) 



for d > 2 
for d = 2 



for d < 2. 



(147) 



Finally, from (|146| ) and ([L47D follows 



lim . 

c(t,t) ] 2 ^ , 



for d > 2 

for d = 2 (148) 
for d < 2. 



These results expose the basic mechanism responsible of 
the behavior of the response function. When looking 
in x space all the k modes are mixed together and the 
existence of one of them which does not equilibrate is 
*q hidden by the density of states as long as d > 2, but 
cannot be canceled for d < 2 and becomes more evident 



As a simple illustration, let us briefly consider the case 
of free diffusion 



which demonstrates quite well the existence of non equi- 
librating degrees of freedom whose visibility depends on 
space dimensionality. In Fourier space ( |l43| ) takes the 
form 



which shows that all the modes with k ^= equilibrate 
while the k = mode executes Brownian motion and 
therefore never equilibrates (23). The linear response 
function is given by R(k,t,t') = exp[— k 2 (t — t')). Inte- 
grating over k and over time, for the integrated response 
function 



the lower is the dimensionality. In particular, ( 148 ) shows 
that the out of equilibrium k = mode does not prevent 
the equilibrium FDT to be asymptotically satisfied for 
d > 2 and also for d — 2, but for d < 2 a deviation from 
equilibrium FDT appears which is increasingly important 
as d — > 1. Interestingly enough, for d = 1 one recovers 
lim t ^oo Tx(t, t w ) = C{t,t)/\f2 as in (§) for the d = 1 
Ising model, recalling that for Ising spins C(t,t) = 1. 
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Appendix A 

In order to compute Xsing(t,t w ) in the d = 1 Ising 
model let us first recall that in the exact solution of the 
model ]2^] the two times and the equal times correlation 
functions are related by 



dt' 



d d k 
(2^* 



-fc 2 (t-t') 



(145) C(\i~j\,t,t') = J2C(\j-l\,t',t')F t _ l (t-t'). (149) 



we find for large time the same pattern of behavior as 
in JT37|) 



,-# 2A° 



d-2 



(4tt 
(47r)~§ log 



(4-r^(^)' 



to 



for d > 2 
for d = 2 

for d < 2. 



where F^ m {t-t') = e (' t ">I l ^ m [^(t-t')} and I n (x) are 
the Bessel functions of imaginary argument. From this 
follows 



% C (\i- j \,t,t') + ^-C(\i-j\,t,t') = 



dt 



(146) 



dt 

E * dT' 



(150) 



The linear response function 
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8 < Sj(t) > h 



(151) 



h=0 



can be cast in the form ll] 



«ij(M') = yfH(<-*'){»(4)) ( 152 ) 

where w(s'j) is the unperturbed transition rate. Rewrit- 
ing the right hand side as 



2 

T ' 



(w(s'j) > Fi-j = ^ | ^(Sj-s^[w;(Sj-) + w(s^)])Fi_ r 



using the parity of Bessel function F x (z) = F^ x (z) and 
the result of Ref. pi 



J2p(u't'\"t) Si = J2s'lFi-l(t-t') 



(160) 



we find 

i m 

= 2 E E ^o^/V, o/Kw, « i o- (i6i) 



- J2 (s'iS'mMs'j) + ™(s'm)])Fi-m } I I ^ I 

and using Glauber evolution equation for C(| i— j |, t', i') 

i-m \,f,f) _ 
dt' 

-2« s m [w( s i) + w(sm)l) form^i , 154 . 
for m = i 



from ( |152| ) and ( |153[ ) one obtains 

TR l^ dC(\j-mlt',t') 
TlUjfr = 2L ^ 

m 

+ B i;i (t,t') 

with 



F%—ffi (t t ) 

(155) 



Bi,j(t,t') = ^2(s' ]S ' m [wis'j) + w(s' m )])Fi- m (t - t'). 
m 

(156) 



Next, inserting (|150| ) in (155) we obtain 



TRi,j(t,t') 



£pC(\i-j\ t tJ) 



g- t C(\i-j\,t,t') 



+ B itj (t,t') (157) 



taking i = j and summing over z this gives 
2 > 



rE^(M') = ^E 



7 C(i,t,*') 



9 



+ B(t,t r ) (158) 



where B(t, f ) = £. B iti (t, t') and C(i, t, t') = C{\ i-j |= 
0, t') is the autocorrelation function. In the general case 
of absence of space translation invariance this quantity 
depends on the site i. Furthermore, from (|156]) 



B(t, t') = J2(«n Ms'i) + w(s' m )])F^ m (t - t') 

■t.rn 

= ^2{Si w ( s i)^2 s 'm) F i- m (t-t') 

i m 

+ 5>>(0 E ^Fi- m (t - t') (159) 



Since the conditional probability obeys the master equa- 
tion El 



^p( W , t I w', = - E S ™ E ^MOpC". ^ I 

(162) 



this gives 



8 



B(t,t') = -J2]f t C(i;t,t' 



(163) 



and putting this result in (158) we finally obtain 

dt 



(164) 



Up to this point the results we have obtained are fully 
general. Let us now specialize to the case of the initial 
condition with a single interface, e.g. oj(t = 0) = [s, = 1 
for i < 0, Sj = —1 for i > 0]. Furthermore, if we take J 
x, also the evolving configuration will contain a single 
interface, namely u)(t) = [si = 1 for i < l(t), s,; = —1 for 
i > l(t)] where l(t) is the position of the interface at the 
time t. If we consider two configurations at the times t, t' 
we have 



J2si(t)si(t')=N-2\l(t)-l(t') 



(165) 



where N is the total number of spins. Taking the thermal 
average 



^2C(i,t,t')=N-2x(t-t') 



(166) 



where x(t — t') = (| l(t) — l(t') |) is the average of the 
absolute value of the displacement of the interface. Since 
this quantity is TTI we may write C(t—t') = C(i, t, t') 
and inserting in (164) we find 



TK{t -t') = 



dC{t-t') 
dt 1 



(167) 
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where K(t - t') = £\ i? M (i,i')- 

Defining X s ing = Q./N) f* w H(t- t')dt' we get 

NT Xsi n g (t - t w ) = [C(t = t w ) - C(t - t w )} = 2 X {t - t w ) 

(168) 

which yields (|94| ) keeping into account that x(t — t w ) ~ 
Appendix B 

Taking for the initial expectations of the auxiliary field 



u(k,0) = 



(169) 



{ u[ki, 0)u(k 2 , 0) = (2<n) d A5(ki + k 2 , 
and using ( |126| ) the unperturbed averages are given by 

d d k 



J(i) = (Vn) 2 
A 



A 



W(fc,t,0)e-A* 



^ 2 (i,o)y (2tt) 



(170) 



and 



u(a?, t)u(x, t') 
A 



= A 



d d k 



-[R(k,t,0)R(k, t',0) 



d d k 



Y(t,0)Y(t',0) J (2ir) d 
where A is the momentum cutoff. Next, using 
d d k 



(2tt) 



and 

we have 
I(t) 



d d k 
(2tt) 



-2k 2 



-2k x 



2Trd(8irx) ~ 



2itdA , , d+2 



t + 



2A 2 



(171) 



(172) 



(173) 



(174) 



From ( p5D 

(ir 2 (f,0) _ 
eft ~ 



2#(t)y 2 (t,0) = -2r e ,F 2 (i,0) + 



u(x, t)u(x, t') 



a 2( 4 ^)rf/2 (t + f + ^) d/2 



(177) 



Using (129) the aging contribution to the correlation 
function is given by 



C ag (t,t') = a(x,t)a(x,t') = m 



2 ll(x, t)u(x, t') 



(178) 



and inserting (131) and (177) this gives 



C ag (t,t') =m\ 



(1 + 1) + 



d/4 



(1 + la 

V2 ~ 2t 



t!_\d/2 
21 ) 



(179) 



which reduces to (134) for to/t — ► 0. 



T /WF 2 (t,0)r 14 
T 

(175) 



and, neglecting the time derivative on the left hand side, 
for large time we find 



Y 2 {t,0) =a 2 [t + 



1 

2A^ 



(176) 



with a = 



3Ad 



4r e „m^(87r)2 



Inserting in (171) we have 



[1] J.P.Bouchaud, L.F.Cugliandolo, J.Kurchan and 
M.Mezard, in Spin Glasses and Random Fields edited 
by A.P.Young (World Scientific, Singapore, 1997). 

[2] G.F.Mazenko, O.T.Valls and M.Zannetti, Phys.Rev. B 
38, (1988) 520. 

[3] L.Berthier, J.L.Barrat and J.Kurchan, Eur.Phys.J. B 11, 
(1999) 635. 

[4] S.Franz and M.A.Virasoro, J.Phys. A 33, (2000) 891. 
[5] L.F.Cugliandolo and J.Kurchan, Phys.Rev. Lett. 71, 

(1993) 173; Philos.Mag. 71, 501 (1995); J.Phys. A 27, 

(1994) 5749. 

[6] S.Franz, M.Mezard, 
G.Parisi and L.Peliti, Phys.Rev.Lett. 81, 1758 (1998); 
J.Stat.Phys. 97, 459 (1998). 

[7] M.Mezard, G.Parisi and M.A.Virasoro, Spin Glass The- 
ory and Beyond (World Scientific, Singapore, 1987). 

[8] A.Barrat, Phys.Rev. E 57, (1998) 3629. 

[9] G.Parisi, 

F.Ricci-Terzenghi and J.J.Ruiz-Lorenzo, Eur.Phys.J. B 

11, (1999) 317. 
[10] C.Godreche and J.M.Luck, J.Phys. A 33, (2000) 1151. 
[11] E.Lippiello and M.Zannetti, Phys.Rev. E 61, (2000) 

3369. 

[12] A. J. Bray, J.Phys. A 22, (1989) L67. 
[13] A.Prados, J.J.Brey and B.Sanchez-Rey, Europhys.Lett. 
40, (1997) 13. 

A first partial account of our results has been published 
in F.Corberi, E.Lippiello and M.Zannetti, Phys.Rev. E 
63, (2001) 061506. 
[15] R.G. Palmer, Adv.Phys. 31, (1982) 669. 
[16] J.Kurchan and L.Laloux, J.Phys. A 29, (1996) 1929; 
C.M.Newman and D.L. Stein, J.Stat.Phys. 94, (1999) 
709. 

[17] For a review see A. J. Bray, Adv.Phys. 43, (1994) 357. 
[18] H.Furukawa, J.Stat.Soc.Jpn. 58, (1989) 216; Phys.Rev. 

B 40, (1989) 2341. 
[19] T.Ohta, D.Jasnow and K.Kawasaki, Phys.Rev.Lett. 49, 

(1982) 1223. 



19 



[20] K.Kawasaki, M.C.Yalabik and J.D.Gunton, Phys.Rev. A 
17, (1978) 455. 

[21] T.Nagai and K.Kawasaki, Physica A 134, (1986) 483; 
A.D.Rutenberg and A. J. Bray, Phys.Rev. E 50, (1994) 
1900. 

[22] R.J.Glauber, J.Math.Phys. 4, (1963) 294. 
[23] L.F. Cugliandolo, J. Kurchan and G. Parisi, J.Phys.I 
France 4, (1994) 1641. 



20 



